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Abstract 

We present a self-consistent calculation of the interlayer asymmetry in bilayer 
graphene caused by an applied electric field in magnetic fields. We show how this 
asymmetry influences the Landau level spectrum in bilayer graphene and the 
observable inter-Landau level transitions when they are studied as a function of 
high magnetic field at fixed filling factor as measured experimentally in Ref. [l[ . 
We also analyze the magneto-optical spectra of bilayer flakes in the photon- 
energy range corresponding to transitions between degenerate and split bands 
of bilayers. 
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1. Introduction 

Bilayer graphene @, H, 0, H, 0, H, Q is one of several graphite allotropes [l(| 
that display unique physical properties determined by the hexagonal symmetry 
of their two-dimensional crystalline structure. The electronic band structure 
of bilayer graphene close to the Fermi energy consists of two degenerate bands 
touching each other at the charge neutrality point (the position of the Fermi 
energy in a neutral system) and two bands split by the interlayer coupling 0, 0, 
[B[ . In bilayer-based field-effect transistors a band gap in the electronic structure 
can be opened using a transverse electric field that breaks the symmetry between 
the layers 0,0,[8, 9], which mayplay an important role in forming the transport 
characteristics of such devices [6(. 

In this paper we analyze self-consistently the interlayer asymmetry parame- 
ter for bilayer graphene field-effect transistors, using the tight-binding approxi- 
mation. The focus of this study is the manifestation of interlayer asymmetry in 
the magnetospectroscopy of bilayers, and we perform a self-consistent analysis 
for bilayer flakes subjected to quantizing magnetic fields, taking into account the 
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Figure 1: (a) Schematic view of the bilayer graphene lattice containing four atoms in the unit 
cell: Al (white circles) and Bl (black) in the bottom layer and A2 (white) and B2 (black) 
in the top one. Also shown are the hopping parameters 70 and 71 used in the tight-binding 
model and interlayer asymmetry u. (b) First Brillouin zone of bilayer graphene with two 
inequivalent K points: K+ and K-. Note the difference between two momentum frames of 
reference: k = + p where k = (k x , k y ) and p = (p x ,Py)- 



possibility that there is finite asymmetry in a neutral structure. When the ma- 
terial is kept at constant filling factor upon the variation of magnetic field (the 
measurement scheme employed in recent experiment [l|), the necessity to vary 
the charge density on the layers causes significant asymmetry and the Landau 
level (LL) spectrum in a strong magnetic field is altered considerably. Our cal- 
culation generalizes the self-consistent analysis Q developed for bilayers at zero 
magnetic field, and the analysis presented here improves the rigour of Ref. [Til ] 
where the variation of interlayer asymmetry on density and its influence on 
Landau level transition energies was estimated by neglecting screening effects. 
Also, we calculate magneto-optical spectra of bilayers in the infra-red spectral 
range that has recently become accessible in optical experiments using cyclotron 
irradiation sources 13, ID, 14 1, and where transitions between degenerate and 
split bands of the bilayer may occur. 



2. The Landau level spectrum in charged bilayers: self-consistent 
analysis 

A schematic view of bilayer graphene (marked with the hopping integrals 
considered throughout this paper) and the Brillouin zone of bilayer graphene 
are shown in Fig. 1(a) and 1(b), respectively. Bilayer graphene consists of two 
coupled hexagonal lattices with inequivalent sites Al, Bl and A2, B2 in the 
first and second graphene sheets, respectively, which are arranged according to 
Bernal (A2-B1) stacking 0. The hexagonal Brillouin zone has two inequiva- 
lent degeneracy points = (£§f,0) (where £ 6 {+, — } and a is the lattice 
constant). Here, we take into account only the nearest-neighbour in-plane and 
interlayer (A2-B1) coupling and the lowest-order terms in the electron band 
energy expansion in electron momentum (determined as the deviation of the 
electron wave number from the corner K + or K- of the hexagonal Brillouin 
zone that is identified below as the centre of the valley), which corresponds to 
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the Hamiltonian: 
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It is written in a basis of sublattice Bloch states W^. — (ipAi,'4 ; B2,'4'A2,'4'Bi) T 
in valley K + and , 3>_ = {i>B2^A\^B\^A2) T in valley and v is related 
to the nearest-neighbour hopping parameter 70 (v = 37^70) an d 7r = p x + ipy 
The interlayer-asymmetry parameter u, which will be at the heart of this study, 
describes the on-site energy difference between different layers in the bilayer. 

In the presence of a magnetic field B perpendicular to the flake, the bilayer 
spectrum splits into Landau levels (LL) . The LL spectrum can be obtained from 
the Hamiltonian in Eq. Jl} using the Landau gauge A = (0, —Bx) and the fact 
that operators n = p x +ip y ~eBx and ir' — p x — ip y + eBx coincide with lowering 
(raising) operators in the basis of Landau functions ip m (x,y) = e lqyy / h (j) m (x + 
q v \g) (where 4> m {x) are the wave functions of a quantum harmonic oscillator), 

h n — , , , + , h 



Kipm = —iT—V2rnipm-i, 7T0O = 0, and i^ijj m = i — \/2(m + l)i/j m +i, (2) 

A_B *B 

where Xb = \J h/eB stands for magnetic length. 

For a symmetric bilayer, u = 0, energies e m of the Landau levels are described 

@, EH EE 13: 



by 
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= i 7? + T(2m - 1) ± v/ 7 4 + 2r 7l 2 (2m — I ) + F 2 ) , m > 1; (3) 



where T = 2^§- and (3 denotes conduction band (+) and valence band (— ) LLs. 

Indices c and s stand for low and high-energy (split) bands and correspond to 
minus and plus signs in front of the square root, respectively. 

Nonzero asymmetry u, caused by a possible difference in electric potential 
energy between the layers, modifies the LL spectrum 0, HH, 0, 21, |2^|. To 



model this effect, we employ a self-consistent theory of the charging of bilayer 
graphene. In order to reproduce the conditions of recent experiments [l[ where 
the transition energy between low-energy LLs was measured as a function of 
a varying magnetic field but for a fixed filling factor, v, we extend the self- 
consistent analysis of Ref. [7| from the zero-magnetic field regime into the 
regime of quantizing magnetic fields, taking into account the possibility that 
there is finite asymmetry in a neutral structure [see Eq. ([4]) below] . The analysis 
presented here improves the rigour of Ref. [l0 | where the variation of interlayer 
asymmetry on density and its influence on Landau level transition energies was 
estimated by neglecting screening effects. 

In particular, we consider a gated bilayer with interlayer separation Co In 
external magnetic field B, a total excess density, n = v^r-, must be induced 
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using the gate in order to keep the filling factor v fixed while changing B. The 
density n is shared between the two layers: n = ri\ + ni where n\ (n 2 ) is the 
excess density on the layer closest to (furthest from) the gate. The difference in 
electric potential between the layers is related to an incomplete screening of the 
gate electric field by the charge en\ on the first layer alone and can be related 
to the unscreened density n 2l 

, m . e 2 c n 2 (is,B) 
u{v,B) = w-\ . (4) 

Here e r is the effective dielectric constant determined by the SiC>2 substrate, 
and w takes into account finite asymmetry of a neutral structure (internal elec- 
tric field due to, for example, initial non-intentional doping of the flake by 
deposits/adsorbates). In our numerical calculations we use e r — 2. 

On the one hand, u influences the LL spectrum via the Hamiltonian in 
Eq. (JT)). On the other hand, its value depends on the charge density n 2 which 
can only be obtained with a full knowledge of the LL spectrum and the wave 
functions corresponding to each LL. Therefore, a calculation of u requires a 
self-consistent numerical analysis. This calculation consists of the following 
steps: for each given B, 5 < B < 20T, and v we choose a starting u, and 
diagonalize Hamiltonian Eq. |T]) to find the LL spectrum and the eigenstates 
with to < M max where M max ~ 300. Then, we sum over all filled Landau 
levels and determine the excess electron densities on each layer. Note that, 



as a nonzero value of u splits the valley degeneracy of the LLs [22j, care has 
to be taken when comparing densities in specific LLs in biased and neutral 
structures, not to confuse levels in different valleys. Finally, using Eq. (|4|) 
we find the asymmetry parameter and, then, iterate the numerical procedure 
to obtain the self-consistent value of u [23j. Note that, for a sufficiently large 
cutoff M max ~ 300, the results were independent of M max . 

The self-consistently calculated values of u obtained for several values of 
the filling factor v are shown in Fig. 2(a) and 2(b) for the case of w = 
and a nominal w — — lOOmeV, respectively. In the case when w = 0, the 
induced intcrlayer asymmetry is antisymmetric with respect to the change of 
the filling factor from positive to negative. This is because changing the filling 
factor from +v to —v corresponds to reversing the applied electric field and 
inducing excess densities — n,— ni and — n 2 and thus reversing the sign of u. 
Also, with decreasing B all curves tend towards u — and u ps — 60meV in Fig. 
2(a) and 2(b), respectively. These values are the results of the self-consistent 
calculation with corresponding values of w in the absence of a magnetic field Q ■ 
Examples of the low-energy LL spectrum for v = +4, w = and for v = +12, 
w = — lOOmeV are shown in Fig. 2(c) and 2(d). To list the LLs in Fig. 2, we use 
three symbols: sm£, where s attributes the LL to the conduction (+) or valence 
(-) band, m is the LL index and £ € (+, — ) identifies the valley (K + or K_) that 
the level belongs to, respectively. The Landau levels m = 0, 1 have no s index, 
as those levels are degenerate when u — 0, H3] ■ The sign of the valley splitting 
of the level sm depends on the sign of u: for u > 0, level sm+ has higher energy 
than level sm— whereas the opposite is true for u < 0. Levels m = 0, 1 behave 
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Figure 2: Top row: results of a self-consistent calculation of the interlayer asymmetry u for 
(a) w = and (b) w = — lOOmeV. Bottom row: the LL spectrum as a function of applied 
magnetic field B for constant filling factor and excess density-induced interlayer asymmetry 
u: (c) v = +4; w = 0, and (d) v = +12; w = — lOOmeV. Solid and dashed lines denote levels 
belonging to K+ and K— , respectively. The line corresponding to the last filled Landau level 
is shown in bold. In these calculations we used v = 10 6 m/s and 71 = 0.4eV. 



differently: in this case the energy e m+ (_) < (>)0 if u > and e m +(-) > (<)0 
if u < 0. The size of the valley splitting of the low-energy LLs increases with u 
and B and for |u| « O.leV, B » 20T [filling factors v = +12, +16 in Fig. 1(a) 
and v = —8, —12, —16 in Fig. 1(b)] is of the order of lOmeV. 

3. Low-energy inter-Landau level transitions and bilayer signature in 
the FIR absorption 

Using spectra similar to those shown in Fig. 2(c) and 2(d), we find the 
energy of the low-energy inter-Landau level transitions for several filling factors 
which mimics the experimental conditions of Ref. [l| (the tight-binding approach 
to this problem has also been adopted in Ref. [ll| where the dependence of 
interlayer asymmetry on density n and its influence on Landau level transition 
energies was estimated by neglecting screening effects). In that experiment, 
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Figure 3: Energy of low-energy inter-LL excitations as a function of magnetic field for w = 0. 
The broken lines are the contributions of individual valleys to the transition energy: black 
dot-dashed and dashed lines denote the transition energy for positive (negative) v at K+ (K— ) 
and K— (-&"+), respectively. The solid black lines show the contribution of both valleys to 
the transition energy, calculated according to Eq. J5) (in this case the transition energy is the 
same for both positive and negative v), whereas solid grey lines depict the transition energy 
in a neutral (u = 0) structure. Note that for v = 8, 12, 16 all black lines are very close to each 
other and difficult to resolve. 



infrared light of energy hw was incident on the graphene bilayer in a strong 
external magnetic field and with a constant filling factor in order to excite charge 
carriers between a prescribed pair of LLs and to probe the energy dispersion. 
According to the selection rules determined in Ref. (and later extended for 
u 7^ [22j). only transitions which change the LL index m by one are allowed. 
Also, as photons provide a very small momentum transfer, we only consider 
transitions between levels at the same K point. Thus, the corresponding low- 
energy transitions for filling factors v = +4, +8, +12, +16 are 1£ — > +2£, +2£ — > 
+3£, +3£ -> +4£, +4£ -> +5£ and for v = -4, -8,-12,-16 they are -2£ -> If, 
— 3£ — ► — 2£, — 4£ — > — 3£, — 5£ — ► — 4£, respectively. However, as transitions 
between the same levels at different K points differ too little in energy to have 
been resolved separately in the abovementioned experiment (in fact, they can 
be only be clearly distinguished in Fig. 3 for the case v — 4), we obtain a single 
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Figure 4: Energy of low-energy inter-LL excitations as a function of magnetic field for w = 
— lOOmeV. Black solid and dashed lines denote the transition energy for positive and negative 
filling factor, respectively. Grey solid lines depict the transition energy in a neutral (u = 0) 
structure. 



transition energy t\ rans for a given filling factor v by comparing the relative 
intensities of the corresponding transition at each K point: 

„ = el ans (K+)I' y (K+) + eUnsiK-ViK-) 

ttrans I" (K+) + I v (K -) [ > 

where e^ rans (K^) and I V {K^) are the transition energy at the iQ point and its 
intensity, respectively. The results obtained for sets of parameters as in Fig. 2(a) 
and (b) are shown in Fig. 3 (w = 0) and Fig. 4 (w = — lOOmeV). We shall first 
discuss the case w = presented in Fig. 3. In this case, for a specified value of B, 
the asymmetry u changes sign with a change of the sign of the filling factor [Fig. 
2(a)], the Landau level spectrum for v and — v are the same but the K points 
have to be exchanged. Therefore, e^ rans (K + ) = e^ ans (K -) , clearly seen in all 
four graphs. Moreover, both transitions have the same intensity and contribute 
equally to e^ rans (black solid line). Comparison with the transition energy for 
a neutral bilayer (grey line) shows that non-zero u decreases the energy of the 
transition. The greater |u| and B, the bigger the difference between excitation 
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energy in symmetric and biased bilayers. However, this difference decreases with 
an increase of filling factor. 

Introduction of parameter w breaks symmetry between the conduction and 
valence band LLs as presented in Fig. 4 for the case of w — — lOOmeV. The 
valence band excitation has greater energy than the conduction band excitation 
at filling factor v = ±4. However, this situation is reversed for higher filling 
factors v — ±8, ±12, ±16 (this reversal was not observed in the experiment 
[1]). For this specific case, w = — lOOmeV, the asymmetry introduced between 
excitations for filling factors v and —v is of the size of 3 — lOmeV. These two 
effects, the reduction of the transition energy with the increase of u and the 
breaking of the symmetry between transitions for positive and negative filling 
factor caused by w, partly account for the disagreement between experimental 
findings and Eq. ([3]) obtained from a tight-binding model for neutral bilayers as 
used in Ref. [H to fit the data. Other investigations 24 , 25| show that additional 
corrections may arise from electron-electron interactions. 



4. IR magneto-optics in bilayer graphene 

In this section, we analyse the optical transition spectra corresponding to 
transitions between LLs in split bands of the bilayers. Here, we use the self- 
consistently calculated interlayer-asymmetry parameter u, and the LL energies 
and wave functions, and we compute the infrared optical absorption spectra 



15|,|26j of right ((B) and left-handed (0) circularly polarized light E w = E u I®/q, 
with I® = (l x — i\y)/y/2 and l e = (la + *lj/)/v / 2 for bilayer graphene in a strong 
external magnetic field. The broadening of the Landau levels is modeled using 
a Lorentzian shape with the same full width at half maximum 7 = 60meV for 
all Landau levels. Numerical results for magnetic field B = 14T and filling 
factors v = 8 and v = 16 are shown in Fig. 5. For the case of w = 0, the 
symmetry of the system demands that the intensity of absorption of light with 
a given polarisation for filling factor v and at the K% point is the same as that 
of the light with the inverted polarisation at the point for filling factor —v. 
This, indeed, is the case for graphs in the left column of Fig. 5, where black 
solid and dashed lines show absorption of right-handed (left-handed) and left- 
handed (right-handed) circularly polarised light for the positive (negative) filling 
factor, respectively. Such a symmetry is broken for the case of w = — lOOmeV, 
for which the spectra for positive and negative filling factors are shown in the 
separate panels, where solid and dashed lines refer to right-handed and left- 
handed polarisation, respectively. In particular, the peak visible for some of the 
spectra at the radiation energy around 0.4eV corresponds to electron excitation 
between the low-energy m = LL and one of the two high-energy m = 1 LLs. 
Its position can be used to determine the value of the coupling constant 71 , and a 
small shift in energy of this peak is due to strong magnetic field and asymm etry 
u. Presented curves show a similar shape as those predicted theoretically IH, 2(| 



and observed experimentally [12j, |13| for IR optical absorption in the biased 
bilayer. 
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Figure 5: Magneto-optical absorption spectra for bilayer graphene in strong external magnetic 
field B = f4T and for filling factors v = 8 and v = 16 (left and right columns, respectively) and 
the case of w = (top row) and w = — lOOmeV (middle and bottom row). For the symmetric 
case of iv = 0, solid and dashed lines show absorption of right-handed (left-handed) and 
left-handed (right-handed) circularly polarised light for the positive (negative) filling factor, 
respectively. For the case of w = — lOOmeV, solid and dashed lines represent absorption of 
right and left-handed circularly polarised light, respectively. 



5. Summary 

In this work, we have considered gated bilayer graphene in external magnetic 
field. We have shown that keeping the filling factor constant results in a break- 
ing of the symmetry between the graphene layers due to the induced carrier 
density. We have calculated the interlayer asymmetry u self-consistently and 
demonstrated its influence on the Landau level spectrum. Using these results, 
we discussed both low-energy and high-energy inter-LL excitations and com- 
pared them to recent experiments concerning optical absorption. In particular, 
we achieved some improvement over the standard tight-binding model for neu- 
tral bilayer in the explanation of the cyclotron resonance experiment probing 
low-energy the Landau level spectrum. 
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